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Informationally complete measurements allow the estimation of expectation values of any opera- 
tor on a quantum system, by changing only the data-processing of the measurement outcomes. In 
particular, an informationally complete measurement can be used to perform quantum tomogra- 
phy, namely to estimate the density matrix of the quantum state. The data-processing is generally 
nonunique, and can be optimized according to a given criterion. In this paper we provide the solu- 
tion of the optimization problem which minimizes the variance in the estimation. We then consider 
informationally complete measurements performed over bipartite quantum systems focusing atten- 
tion on universally covariant measurements, and compare their statistical efficiency when performed 
either locally or globally on the two systems. Among global measurements we consider the special 
case of Bell measurements, which allow to estimate the expectation of a restricted class of operators. 
We compare the variance in the three cases: local. Bell, and unrestricted global — and derive condi- 
tions for the operators to be estimated such that one type of measurement is more efficient than the 
other. In particular, we find that for factorized operators and Bell projectors the Bell measurement 
always performs better than the unrestricted global measurement, which in turn outperforms the 
local one. For estimation of the matrix elements of the density operator, the relative performances 
depend on the basis on which the state is represented, and on the matrix element being diagonal or 
off-diagonal, however, with the global unrestricted measurement generally performing better than 
the local one. 

PACS numbers: 



I. INTRODUCTION 

One of the most distinctive features of quantum mechanics is that the information on the state of the quantum 
system cannot be retrieved in a single measurement Ij . Among the different methods for retrieving information on the 
state — including quantum tomography, state discrimination, and quantum state estimation Q — the informationally 
complete measurement ^ |1| (shortly named infocomplete in the following) is the most versatile, allowing one to 
obtain the expectation (O) of any operator O of the system, and with the additional possibility of adapting the 
processing of outcomes, depending on the quantity to be estimated and on the set of states on which the measurement 
is performed. A particularly interesting case is when the measurement is covariant 6, 7] with respect to a group of 
physical transformations. This means that there is an action of a group on the probability space which maps events 
into events in a covariant fashion, namely that when the quantum system is group-transformed, the probability of 
the given event becomes the probability of the transformed event. This situation is very natural, and occurs in many 
applications (e .g. heterodyne detection 0, rneasurement of the spin-direction P|, transmission of reference frames 
[lOj, estimation of unitary transformations etc.) 

In this paper we study covariant infocomplete measurements for bipartite quantum systems. The problem of 
optimization of data-processing will be solved using the method of frame theory 0, |3 , which recently has been 
shown 0, 01 to be perfectly suited to optimize covariant infocomplete measurements. Dealing with bipartite (or 
multipartite) quantum system, the study of the possibility to perform information-theoretic tasks locally on non-local 
states has become a new paradigm in the field of quantum information. In this respect consider for example the study 
of local implementation of non-local quantum gates [1^], telecloning [l7| . remote manipulation and preparation of 
quantum states [T^ , and local discrimination of non-local states [T9l | . The problem of comparing the performances of 
local and entangled estimation schemes has not been addressed yet, and this is precisely the purpose of the present 
paper, where we compare the statistical efficiency of local infocomplete measurements with that of nonlocal ones. In 
particular, the possibility of having a maximally entangled infocomplete measurement — or Bell quorum — has never 
been considered. Here we will show that, although an infocomplete Bell measurement strictly does not exist, however, 
Bell measurements can be used to estimate the expectation of any operator O having both partial traces Tri [O] 
and Tr2[0] proportional to the identity, e. g. maximally entangled operators and traceless factorized observables. 
This possibility is interesting also in view of the relevant problem of the classification of measurements according to 
separability criteria, which still remains unsolved. 

The statistical efficiency of an infocomplete measurement is quantified by the noise in the estimation of the operator 
expectation (O). This generally depends both on the state and on the operator O to be estimated. However, for 
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unbiased estimations the squared expectation is independent of the measurement, whereas if the measurement is also 
fully covariant (i. e. covariant under the full group of unitary transformations) the square of the estimation averaged 
over the ensemble of possible input states depends only on the measurement: this greatly simplifies the evaluation of 
the variances, and allows a simple comparison of statistical efficiencies of different infocomplete measurements. 

The paper is organized as follows. After summarizing some useful notation and elements of the Theory of Frames in 
Sectionini in Sect. Illll we give a general method for optimizing the data-processing. In Sect. II VI we consider covariant 
infocomplete measurements for bipartite quantum systems, in which case the optimal data-processing is provided by 
the so-called canonical dual frame. In Sections lVllVIl and lVIll we derive the minimal noise (averaged over an ensembles 
of covariant states) for factorized, global, and Bell infocomplete covariant measurements, and finally we conclude the 
paper in Section IVlIII bv comparing the resulting statistical efficiencies of the different types of measurements. 



II. NOTATION AND ELEMENTS OF THE THEORY OF FRAMES 



An infocomplete measurement is characterized by a Positive Operator- Valued Measurement (POVM) that spans 
the whole space of linear operators. The Theory of frames [l^ provides simple and powerful tools to establish 
whether a set of operators is complete in the whole operator space, and in addition gives a simple algebraic rule for 
constructing all possible expansion coefficients — representing all possible data-processing of measurement outcomes — 
in terms of the so-called dual frames. Upon denoting with {Pi} the elements of an infocomplete POVM dual 
frames correspond to sets of operators {Di} in terms of which we can write the operator expansion as follows 

N 

^Tr[i?j0]P, = 0, (1) 

1=1 

where N denotes the number of outcomes of the POVM (in the following we will not specify it anymore in the sum 
limits). When considering expansions of operators over a Hilbert space TL it is convenient to exploit the natural 
isomorphism between operators O on TL and vectors | O)) in 7i®^, defined through the equation 

\0)) =Y,{MO\n)\m)\n) . (2) 

We will make repeated use of the following identities po| 

A®B\C))^\ACB^)), (3) 
^:T^[\A)){{B\]^A-B\ (4) 
TT2[\A))iB\]=AB\ (5) 

where r and * denote transposition and complex conjugation with respect to the given fixed basis in Eq. The 
main results from frame theory are the following 13). First, a the POVM {P,} is infocomplete if and only if the frame 
operator 



(6) 



is invertible. A dual frame {Di} is given by the following one (usually referred to as canonical) 

|A))=F-^|P,:)). (7) 
In general there exist infinitely many alternate duals {D[} , which can be obtained from the canonical one as follows 

D[=D, + Y,~Y.^v[P}d,]Y,, (8) 

3 

for arbitrary set {!}}. Finally, whenever a POVM is not complete, the span of its elements coincides with the 
support (the orthogonal complement of the kernel) of the operator F. Then the operators in such a subspace can 
be reconstructed through Eq. (^), where the canonical dual is now defined through the Moore-Penrose generalized 
inverse F^ of F , which corresponds to invert F on its support. 
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III. OPTIMIZATION OF THE DATA-PROCESSING 



In the following wc will restrict our attention only to finite dimension d < oo of the system Hilbert space. Apart 
from the case of POVM {Pi} with exactly (P linearly independent elements (whence with discrete sample space), 
for infocomplete POVM's there always exist infinitely many alternate duals. This feature provides a wide freedom 
in choosing the dual frame, which can be exploited in order to optimize the data-processing, e. g. to minimize the 
variance in the estimation. Therefore, suppose that an experimenter is collecting statistics using an infocomplete 
measuring apparatus — e. g. in a quantum tomographic setup — with the aim of estimating the expectations of any 
desired operator. The noise in the estimation of the operator O can be evaluated by 

50\p-P) = |Tr[i?l0]pTr[P,p] - |Tr[Op]p , (9) 

■i 

which clearly depends on the state p. Notice that we consider generally complex operators O, i. e. we can include the 
case of external products O = of orthonormal basis, which are needed to estimate the matrix elements pij of 

the density operator p. In a Bayesian scheme the noise (O can be minimized for a given prior probability distribution 
d^{p) of the input states p, corresponding to the prior state p = J dp{p)p. The dual frame {Di\ can be chosen in 
order to minimize such a noise, which in turn corresponds to minimize the following norm 

^2 1^1' =^E|Tr[^l0]pTr[P,p], (10) 

i 

where |c) = Yld^i Tr[i?|0]|ei) is a vector in C''^ , and {ci} denotes the canonical basis of C^. We will now state a 
useful condition for optimality. Let us consider the linear mapping A : — > 7i®^ 

A|c)=5]c.|P.)), (11) 

i 

where Ci = (ei|c). The matrix elements representing A in the bases |ei) G and |rn)|n) e 7i®^ are given by 

Ajnn,z — {Pi)mn • (-^^) 

We can define F to be a generalized inverse of A, shortly g-inverse, if A = AFA holds (in the following, for all properties 
of the various kinds of g-inverse see Ref. (23]). We now show that any g-inverse F must have matrix elements of the 
form 

where {Di} is a dual frame for {Pi}. In fact, (FA)y = Tr[Z3|Pj], and consequently the condition for F to be a 
g-inverse can be written as {ATA)mn.j = X)i(-P'i)mnTr[-DjP,]. Then, by informational completeness of the POVM 
{Pj}, the reconstruction formula ^iT^lOjOjPi = O must hold, and necessarily {Di} is a dual. Moreover, since 
{Di} is a dual, than the g-inverse is also reflexive , namely FAF = F, as can be easily checked from the expansion 
(FAF)i_,„„ = Ti [DiPj]{Dj)n7n — {D*)mn = (r)i, mn- In addition, the g-inverse F is also a least square inverse, since 

||AFO — 0|| = 0, as a consequence of the expansion (AFO)„i„ = Tr[I?|0](Pi)„i„ = (0)m„. Therefore, summarizing, 
any dual frame {Di} corresponds to a reflexive and least square g-inverse F of A through Eqs. p2l) and (|13|l . 
Now, the quantity we want to minimize is ||FO||7r, where the norm || • Itt is defined through 

l|c||^ = (cklc) , (14) 

TT being the diagonal matrix (in the canonical basis e^) with entries tth = Tr[Pip] (for finite dimension the POVM is 
trace-class). Notice that when Tr[Pip] does not depend on i (as in the case where the average over the prior distribution 
of states gives p cx / and the POVM is covariant so that Tr[Pip] oc Tv\UivUl] — Tr[i/]) then tt oc I. The quantity to be 
minimized is then simply ||cp = (c|c), and the solution to the minimum- norm problem is provided by any matrix F 
such that FA = A^^F^ = A^F^^FA. Along with the reflexivity and least-square properties, the minimum-norm condition 
uniquely determines the g-inverse in terms of the Moore-Penrose pseudo inverse F = A^. Since (FA)y = Tr[£'JPj], 
the previous identity is equivalent to 



k 



(15) 
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Now it is easy to check that the canonical dual Di, defined as \Di)) = F ^\Pi)), satisfies the previous condition, as 

((A|P,)) = mm) = J2{{Pm)){{Dk\P,)) = mF-'\P,)) . (16) 

k 

In the general case in which tt 9^ /, analogous proof as in the previous case P3| leads to the following condition for 
minimization of ||rO||7r 

ttFA = A^rtTT ^ A^rtTrFA . (17) 

The results of the present and the previous sections also apply in the case of continuous POVM by suitably replacing 
the discrete index i with a continuous one, and sums with integrals with some care about summability of integrals in 
the case of noncompact groups. 



IV. INFOCOMPLETE COVARIANT MEASUREMENTS FOR BIPARTITE SYSTEMS 



In the following we will consider quantum measurements that are covariant under the action of a group G of unitary 
operators Ug, g ^ G, on the Hilbert space 7Y®^ of a bipartite system, H denoting the Hilbert space of the two identical 
systems, for finite dimension d — dim(7i). As well known 6], the POVM of the measurement has the form 



dgPg^dgUg^Ul 



(18) 



where ^ is a suitable positive operator such that the POVM is normalized, and dg denotes a (suitably normalized) 
Haar invariant measure on the group (for what we need in the following, the group is unimodular, which guarantees 
that an invariant measure always exists) . For later convenience, we will normalize the Haar measure over the group as 
Jf^dg = d {d is the dimension of the Hilbert space on which the group is represented) . For infocomplete measurements 
every operator can be expanded over the POVM density, namely Pg spans the whole linear space of linear operators. 
Clearly, all expansions for bounded operators must be square summable over the group G. However, since we are 
considering only finite dimensions, all operators are bounded, and since we consider only compact groups (which then 
admit normalizable invariant Haar measure), all group integrals are bounded too. In the present case, the noise in 
Eq. @ can be rewritten as 



SO^{p;P)= [ dg\TT[DlO]\'TT[Pgp]-\TT[Op]\' 
Jg 



(19) 



We will consider the following classes of pure states with uniform prior probability distribution: a) all pure input 
states; /) all factorized pure input states; e) all maximally entangled input states (also called EPR states). The 
averaged noise in all three cases can be evaluated with the following integrals 



6aO^[P] = i / dhSO^ (Uh\0))mulP 
'^/O't^l = i / d/i / dh' 5O'({Vh®Vh'm{0nv;i®Vl);P 

" JsV(d) JSV{d) ^ 

6eO^[P] = i / dhSO^ ({Vh®I)J{Vfl^I);P) , 

" JsV(d) ^ ' 



(20) 
(21) 
(22) 



where X — ^|/))((/|, while |0)) and |0) arc arbitrary reference pure states in 7i®^ and 7i, respectively. Notice that 
the average of |(0)P = |Tr[Op]p over input states does not depend on the POVM, but only on the set of states over 
which the average is estimated. Upon denoting by overline the average over input states, one has 



m\l = ^ / d/.Tr 



SU(d2) 



®2 



d2 
1 



•m{d) JBV{d) 



\{0)\} = / dh dh'Ti {Vh^Vh'f^\0){0\''\v;i®V,l,)''\O^O^) 



\{0)\l = cl/iTr {Vh®I)®^T®\vI®I)®\0®0^) 

d Jsv{d) L 



(23) 
(24) 
(25) 
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These integrals can be evaluated by exploiting the following identities (which are corollaries of Schur's lemmas) 

/ dgUgXUl=Tv[X]I , (26) 

/ dff UfXUl^' = -^Tr [P^X] PP + -i-Tr [P^X] P^ , (27) 

where Pg' and P^ denote the projection on the symmetric and antisymmetric subspaces of /C*^^, respectively. The 
result is 

4 



l(0)l? = ^2(^^Tr[(P|^,3®P|^,J(0®Ot)] ^ 

^ ^ d3(d+l) ^^' ^(^^""13 ® 24)(0 ® O^)] , (28) 

where denotes an operator acting on Hi <S) Hj . Finally, recalling that Ps = ^{I + E) where E is the swap operator 
E\(j))\'ip) = \ip)\<j)), and that Tr[E{A ® B)] = Tr[AB], we get 



1 




M) 


1 




d2(d + 


1)2 


1 




2d^{d- 


f 1) 



(29) 



= 4,1^2 (Tr [m + |Tr[0]|2 + Tr [|Tri[0]|2 + |Tr2[0]p]) , (30) 



= ^5:«77-TT + |Tr[0]p + Tr [\Tt,[0]\' + |Tr2[0]p]) , (31) 



Notice that |(0)|2 = ±ti|(0)|2. 

One can easily show that the first integral in Eq. (|19|) is independent of the input state ensemble, so that the first 
term of the noise depends only on the POVM, whereas the second term depends only on the input ensemble. This 
can be verified using identities and (|?7|l . One has 

SOl[P] Jjg\Tr[DlO]\'Tr[P,]^](m, (32) 

where x = a, f, e, and a dual that optimizes one of these noise parameters optimizes all of them. According to Sect. 
mil we are in the situation where the canonical dual is optimal, namely it provides the optimal processing function 
minimizing the noise (|19() . In the following section we will then consider processing functions obtained from the 
canonical dual. 



V. PRODUCT OF LOCAL INFOCOMPLETE MEASUREMENTS 

As the first example of infocomplete POVM we consider 

Pl°^ = Ugi^Ul ® U,yul , (33) 

where the elements g,h belong to SlLJ((i), and ly,!/' are pure states in C^. Such a POVM describes a measurement 
that can be performed locally by two separate parties, with classical communication needed in order to evaluate the 
processing function. 

It can be easily shown by Schur's lemma that Pg ® Ph is actually a POVM. The canonical dual can be written as 

M 

Dg,h = [{d + l)UgvUl -I]® [{d + Iphiy'ul - I] . (34) 
The noise 6xO^ in the evaluation of the expectation value of an operator O is given by 

5.0^[P'°^] = ^[ dg f dh\Tr[Dlp]\^-J{Om- (35) 
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Substituting the expression for the canonical dual and exploiting the identities in Eqs. (|26|l and (|27() . we obtain 



1 



- (d+ l)Tr [|Tri[0]p + \Tv2[0]\'] } - \{0)\l . (36) 

VI. GLOBAL INFOCOMPLETE MEASUREMENT 

We will now consider the POVM 

Pt" = U,,.Ul , (37) 
where g now belongs to §U((i^), is a pure state in ® C^, and /§u(^2) = (P. The canonical dual set is given by 

i^g = (d^ + l);7<,i.{7] - / . (38) 
The average noise over all pure states can be evaluated as in the previous section, and one obtains 

^O2[pgiob] = 1 {(d^ + l)Tr[|On _ |Tr[0]p} - \{d)\l • (39) 

VII. BELL MEASUREMENT 

Finally, we will consider the Bell POVM 

pBeii ^ c/^ ^ J I J)) ((J \ul(g>I = d{Ug {Ul <E> I) . (40) 

with g belonging to SlLJ((i). Using Eq. H26(l it is easy to verify that Eq. H4Q|I actually defines a POVM. On the other 
hand, such a POVM is not infocomplete since operators whose partial trace is not proportional to the identity cannot 
be spanned. This can be seen directly from the reconstruction formula (Q. We can evaluate explicitly the frame 
operator as follows 

F = f dg\P,)){{Pg\ 

= <f I dg Ug®I®U; ® /(Ii2 ® ImWI ®I®Ul®I 

= d|ji3®l24 + ^^J^(/l3-2:i3)®(/24-2'24)| , (41) 

where denotes the projector I acting on Tii ® Hj. Indeed, F has a nontrivial kernel. The Bell POVM in Eq. 14U|) 
can be used, however, to obtain the expectation values of any operator of the form 

O = aoo/ + ® ^3 ' (^^2) 

where Vi and Wj are orthogonal bases for the operators with Vq = Wq = / (and hence with Vi and Wj traceless for 
i,j ^ 0). These are exactly the operators in the support of F. Notice also that any maximally entangled state is of 
the form ll^ . 
Since one has 

^^^1 f/.))l2| C/;)>34 = ^1 C/,)>12| ^;))34 - ^1 /)>13| /))24 , (43) 



the canonical dual is given by 



, 

IF 



Dg = —r- 1 Ug)) {{Ug I - —j^I . (44) 
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As long as operators O of the form H42() are considered, the average noise can still be evaluated through Eq. (|32|l , and 
is given by 

JO^fpBoU] ^ 1 /• dg\TT[D,0]\'-\{dm- (45) 

" JSV(d) 



Evaluation of the integral gives the following result 



1 



^O^I^n = :^{(d'-l)Tr[|Op] + |Tr[0]p} 



d3 



TT[\Tr,[Or + \TT,[0]\']-\{0)\l. (46) 



VIII. COMPARISON AND CONCLUSIONS 



We now compare the three infocomplete POVM's in terms of their statistical efficiency, namely we compare their 
variance H19|l . We have 



<50^[pioc]„^02[pgiob] ^ 2^Tr[|On + i|Tr[0]n-^Tr[|Tri[0]p + |Tr2[0]n , 



da 

(502[pi°c] _5o2[pBcn] = ^(ltl)Tr[|0|2]-^Tr[|Tri[0]|2 + |Tr2[0]p] . (47) 

Clearly, for operators O having Tri[0] = Tr2[0] = (e. g. for operators of the form O = A B with both A and B 
traceless) one has 

For off-diagonal matrix elements in a factorized basis, namely with O ~ 'S> \n){m\ and i ^ j and/or n ^ 

the inequalities (|48|l also hold true, and, specifically, for diagonal matrix elements one has J^O^i^'^'"''] = '5^0j;[P'°'^]. 
Also for maximally entangled projectors O = the variances are bounded as in Eq. H48|) . whereas for 

O = i| V)){{W I and Tr[H/tp-] = q one has 

and generally the expectation of O cannot be estimated using the Bell POVM, unless W'^V cx /. This means that 
for matrix elements in an orthogonal Bell basis | Vnf) where 14, are traceless (apart from the identity) and form a 
also a (projective) representation of a group (such as the usual orthonormal basis with Pauli matrices for d = 2, or 
shift-and- multiply operators in dimension d > 2), the Bell POVM is optimal for diagonal matrix elements, whereas 
for off-diagonal matrix elements it cannot be used, and for these it is better to use the global measurement than the 
local one. 

We conclude that for factorized traceless operators, for Bell projectors, and for the density matrix estimation in 
factorized basis the Bell measurement always performs better than the unrestricted global measurement, which in 
turn outperforms the local one. For the diagonal matrix elements in the factorized basis the Bell measurement cannot 
be used, and global and local POVM's perform equally well. Finally, for matrix elements in an orthogonal Bell basis 
the Bell POVM is optimal for diagonal matrix elements, whereas for off-diagonal matrix elements it cannot be used, 
and for these elements the global measurement performs better than the local one. 

A further interesting study continuing the present analysis is the evaluation of the noise for infocomplete measure- 
ments with minimal number of outcomes, comparing local versus global measurements. 
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